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I. INTRODUCTION 

The foundations of quantum theory have been known for 
several decades, but the relation to classical physics is still a 
topic of current research. In quantum optics, the notion of 
nonclassicality caused long-lasting discussions. The quan- 
tum state of a radiation field is often examined by means 
of photodetectors which measure normally ordered field cor- 
relation functions. The latter are properly described by 
the quasiprobability distribution or P function of Sudarshan 
and Glauber fl] |2). Following Titulaer and Glauber J5], 
"states with positive P functions ... are ... possessing clas- 
sical analogs." The other way around, a quantum state is non- 
classical if the P function does not exhibit the properties of a 
classical probability distribution, cf. e.g. 

Any quantum state of a harmonic oscillator can be given as 
a quasimixture of coherent states \a), 

p = J d 2 aP(a)\a){a\, (1) 

where P(a) is the P function mentioned earlier, cf. (T||2). 
The coherent state is known to be that quantum state which 
is most closely related to the classical behavior of an oscilla- 
tor. Its P function is formally equivalent to the deterministic 
classical phase-space distribution, representing a single point 
in phase space. If the P function has the properties of a clas- 
sical probability density, P(a) = P c \(a), the state is a true 
classical mixture of coherent states. Hillery has shown that 
the coherent states are the only pure quantum states having a 
non-negative P function Q. Hence, for any classical mixture 
of coherent states the P function exactly reflects the classi- 
cal behavior of the oscillator in phase space - including its 
free evolution. The failure of the interpretation of P(a) as 
a probability density, P(a) ^ P c \(a), is intimately related 
to the quantum superposition principle, thus it most naturally 
displays the quantumness of any quantum state. 

However, in general the P function can only be understood 
as a generalized function which is often not accessible. For 
this reason, different representations of a quantum state are 
considered. An often used one is the Wigner function |6|, 
which also covers the full information on the quantum state. 
A generalization yields the set of s-parameterized quasiprob- 
ability distributions [7]. By fixing the parameter s, different 
quasiprobabilities are obtained. If one of these functions vi- 
olates the requirements of a classical probability distribution, 



the given state is nonclassical. Unfortunately, this set of func- 
tions does not reveal all nonclassical effects in terms of regular 
functions: For a squeezed state, they are either non-negative or 
highly singular. In order to develop quasiprobabilities to un- 
cover nonclassicality in general, the generalized quasiproba- 
bilities of Agarwal and Wolf will be a powerful foundation JS] . 

Another general representation of a quantum state is its 
characteristic function, defined as the Fourier transform of a 
given quasiprobability. Its advantage lies in the fact that it is 
always a regular function, even the characteristic function of 
an irregular nonclassical P function. Useful nonclassicality 
conditions have been derived (51 [TO) and applied in experi- 
ments lfTTHT3l . However, for a full characterization of non- 
classicality one needs to check an infinite hierarchy of condi- 
tions, which may be a cumbersome procedure. We will use 
them as the starting point for our examination. 

In this article we introduce regular quasiprobabilities with 
the aim to uncover all types of nonclassical effects by their 
negativities. A distribution of this type, to be called nonclassi- 
cality quasiprobability, belongs to the set of the Agarwal- Wolf 
quasiprobabilities. For our purposes the filter functions occur- 
ring in the latter must obey specific constraints. We study the 
properties, which are needed to make the filters useful for ex- 
perimental applications, and we show how to construct them. 

The article is structured as follows. In Sec. II we introduce 
the requirements for general nonclassicality filters and discuss 
the relation to previously known filter procedures. Section III 
is devoted to nonclassicality quasiprobabilities, which, in ad- 
dition, contain full information about the quantum state. The 
method is illustrated for the example of a squeezed vacuum 
state. In Sec. IV we briefly summarize our results. 



H. NONCLASSICALITY FILTERS 

A. Characteristic functions and Bochner's theorem 

Let us now consider the possibility of getting general in- 
sight into the properties of the P function in an experiment. 
The characteristic function <&(/?), defined as the Fourier trans- 
form of P(a), can be sampled by balanced homodyne detec- 
tion, cf. lfT3l[l4ll . From a set of quadrature data {xj [v]}|Li at 
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some fixed phases ip, it can be estimated by 
i N 

= J_ V" e i|^k 3 [7r/2-arg(/3)] e |/3| 2 /2^ ( 2 ) 

Thus we have direct experimental access to the characteristic 
function (f>. It may be a rising function of whose Fourier 
transform only exists as a highly singular distribution lfl3l . 
The standard deviation of $ is, for a given sample of data, 
bounded by lfT4l 

e l/3| 2 /2 

am/3)} < -—. (3) 

Interestingly, also, the characteristic function $(/3) of any 
quantum state is bounded by e'^' I" 1 |[T5l . 

If is not square integrable, its Fourier transform is 
highly singular. In such cases, the nonclassicality of the 
quantum state can be identified via Bochner's theorem lfl6l . 
The function is in general continuous, with $(0) = 1, 
$(—/?) = $*(/?). The P function has the properties of a 
probability density if and only if for all positive integers N 
and complex ft, ...,(3 N , the matrix ($(ft - Pj))ij=\,...,N 
is positive semidefinite. This leads to an infinite hierarchy 
of nonclassicality conditions IfTOl . which in practice cannot 
be examined completely. However, for N — 2 we obtain a 
simple inequality, which is valid for all (3 and necessary for 
classicality 0: 

|*(/8)| < 1- (4) 

The violation of this inequality can be used to experimentally 
demonstrate the nonclassicality of a quantum state lfTTT - IT3ll . 

In cases when the inequality |4} is fulfilled, one cannot di- 
rectly infer classicality, but there is a chance that the charac- 
teristic function is square integrable. A prominent example 
are the photon-added thermal states fP71 . Then one can per- 
form the Fourier transform to obtain the P function and check 
nonclassicality by its original definition lfl4ll . Severe problems 
occur to identify nonclassicality if the characteristic function 
satisfies but is not square integrable. 

B. Filtered P functions 



such that its Fourier transform, Pq([3;w), is a well- 
behaved function. Since $(/?) and its standard devia- 
tion (cf. Eq. ([3])) are bounded by e^l"/ 2 , we need that 
O u ,(/3)e'^'"/ 2 is square integrable for all w. This en- 
sures that the method is universal: it applies to any 
quantum state and to realistic experimental data. 

(b) Non-negativity. To detect nonclassicality of unknown 
quantum states by negativities in the regularized func- 
tion Pq(/3;iu), the latter shall be non-negative for all 
classical states. Equivalently, the filter Q w (/3) shall not 
cause additional negativities in the regularized function 
Pq. This requires that ft w ((3) itself has a non-negative 
Fourier transform. 

(c) Completeness with respect to the nonclassicality of the 
P function. The parameter w represents the width of 
the filter. It may be introduced by a scaling transform, 

Sl v ,(0) = n 1 (J3/w). (6) 

For an infinitely wide filter, the Pq function approaches 
the original P function. This requires, for all (3, that 

lim $ n (M = *03), (7) 

w— >oo 

or equivalently, £l w (0) = 1 and lim. l „_ i . 00 £l w (ft) = 1. 

The most simple example of such a filter is a 
two-dimensional triangular filter, Q, w (f3 r + i(3i) = 
tri(/? r /w)tri(/3i/w), where tri(x) = 1 — \x\ for \x\ < 1 
and tri(x) = elsewhere. Since this function has compact 
support for all iv > 0, it satisfies the condition (a). Further- 
more, it obeys the constraints (b) and (c), since the Fourier 
transform of the triangular function is non-negative and 
lim^^oo fl w (/3) = 1, respectively. 

This example clearly shows that there exist filters which sat- 
isfy all our requirements. Most interestingly, they can be used 
to detect nonclassicality of any nonclassical state. The other 
way around, the negativities are uniquely caused by the non- 
classicality of the state, not by the filter. For all nonclassical 
states, we can find a regularized function Pq which displays 
negativities. We refer to such filters fl w (/3) as nonclassicality 
filters. For the proof of their general properties, we refer the 
readers to theorem 1 in Appendix [A] 



Let us now develop a simple and general method for iden- 
tifying the nonclassicality of a quantum state under realistic 
experimental conditions. It is based on filtering of the charac- 
teristic function, 

$n(/3;w) = $(f3)n w (P), (5) 

with a filter function il w ((3), which we will allow to depend 
on a real parameter w. The filter shall satisfy the following 
specific properties: 

(a) Universality. For any quantum state, the filtered 
characteristic function $si(/3;w) is square integrable, 



C. Relation to known filtering procedures 

Filtering procedures of the P function having the structure 
of Eq. |5]l are already known. However, there is no procedure 
known that fulfills all the requirements (a)-(c). Let us briefly 
consider such filtering approaches together with their short- 
comings for nonclassicality detection. Note that the following 
approaches had not been designed for that purpose. 

By choosing Q s (0) — exp((s — l)|/3| 2 /2) as a family of 
filters, we consider the s-parameterized quasiprobabilities Q. 
For s = 0, we get the Wigner function; for s = 1, the P func- 
tion; and for s = — 1, the Q function. It is obvious that such 



3 



filters do not fulfill the universality condition (a) for s > 0. 
Therefore, they are not capable of regularizing the P function 
of an arbitrary state for s > 0. There exist nonclassical states 
which do not possess a regular s-parameterized quasiproba- 
bility showing negativities. Squeezed states are a prominent 
example. If their nonclassical effects would already be dis- 
played for s — 0, they could be observed as negativities in the 
Wigner function. However, the Wigner function of a squeezed 
state is always non-negative. 

Another filter was considered by Klauder ifTHl . He showed 
that appropriate filtering of the P function may lead to an 
infinitely differentiable regular function. This filtering was 
recently applied to regularize the P function of a squeezed 
state lfT9l . However, since Klauder's filtering does not obey 
the non-negativity condition (b), the corresponding negativi- 
ties of the regularized functions are not uniquely related to the 
nonclassicality of the considered quantum state. 

Last but not least, a very general approach to define 
quasiprobabilities and operator ordering was introduced by 
Agarwal and Wolf JSJ . This may be considered as a general 
background of our considerations. Since the authors' aim was 
to provide general methods, they did not consider constraints 
of the type needed for the nonclassicality filtering. 



III. NONCLASSICALITY QUASIPROBABILITIES 
A. Filters for quasiprobabilities 



Since Vt w {j3) is continuous, the sequence of functions con- 
verges for all (3 point wise to 1 when w — >• oo: 

lim fl w (p) = lim 9,{fi/w) = fi(0) = 1. (11) 

Hence these functions satisfy all criteria for being a nonclas- 
sicality filter. Since 0«,(/3) has no zeros, the regularized func- 
tion Pq contains all the information about the quantum state, 
and consequently, it represents a generalized phase-space 
function in the sense of Agarwal and Wolf [8 |. Therefore, for 
any nonclassical state, one can find a regular quasiprobability 
distribution which displays the nonclassical character by its 
negativities. We refer to such distributions as nonclassicality 
quasiprobabilities . 

The experimental implementation of the procedure to iden- 
tify nonclassicality of a general and unknown quantum state 
is straightforward: 

(A) Sampling. Direct sampling of the function $(/?) from 
experimental data, Eq. d2j, and estimation of its stan- 
dard deviation, cf. Eq. (pT 

(B) Filtering. Choose the set of nonclassicality fil- 
ters S} w (/3), for example the autocorrelation filters in 
Eqs. ^ to ([To}. Multiply the sampled $(/}) with the 
filter of width w: $n(P;w) = $>(0)n w {0). The single 
parameter w is used to optimize the statistical signifi- 
cance of the nonclassical effects to be visualized. 



Filters with compact support can be easily applied to ex- 
perimental data. However, one loses information about the 
quantum state, such that the latter cannot be recovered com- 
pletely from a filtered P function. To overcome this problem, 
one has to use invertible nonclassicality filters. They have to 
meet the criteria of Agarwal and Wolf (8], in particular, hav- 
ing no zeros anywhere, in order to preserve all information 
about the state. 

To our knowledge, no simple examples for such filters are 
known, but they can be constructed in the following way: Let 
us assume that some positive continuous function ui(j3) sat- 
isfies Lj(—j3) = uj{(3) and decays sufficiently fast; that is, 
w(/3)e"'' 3 ' j s square integrable for any u > 0. For example, 
one may choose 



W (/3) = expH/3| 4 ). 
It is easy to see that its autocorrelation function, 



ncs) 



i 

AT 



u;(/3>(/3 + /3V 2 /3', 



(8) 



(9) 



with Af — J |w(/3)| 2 d 2 f3, is positive and satisfies f2(-/3) = 
Moreover, we find that O(0) = 1 and the Fourier trans- 
form of an autocorrelation function is always non-negative. 
Finally, it decays sufficiently fast, so that r2(/3)e u ' /3 l is square 
integrable for any u > 0(cf. lemma 1 in Appendix [B|. 
Now we define a set of functions by 



n w (i3) = n(/3/w), w>o. 



(10) 



(C) Fourier transform. Calculate the Fourier transform Pq 
of $fj and its statistical error. If it displays statistically 
significant negativities, the state is clearly nonclassical. 
The wider the filter, the more nonclassical effects are 
visible in the regularized function Pq, which is only 
limited by the increasing sampling noise. 



B. Example: Squeezed vacuum state 

For illustration, let us consider a squeezed vacuum state, 
described by a characteristic function 

$(/3) = exp{-(/3 + /3*) 2 V x /8 + [fi - /T) 2 V« + l/?| 2 /2}, 

(12) 

where V x and V p are the variances of two orthogonal quadra- 
tures and V x < 1 < V p . The P function is highly singular: 
It is composed of derivatives up to infinite orders of the Dirac 
5 distribution. Hence it is extremely difficult to verify the non- 
classicality of a squeezed vacuum state in this general sense. 

All s-parameterized quasiprobability distributions are ei- 
ther Gaussian or highly singular, and therefore none of them 
has negativities which can be directly reconstructed. Let us 
now consider a squeezed state with V x = 0.2 and V p = 5.0, 
which can be experimentally realized. We construct the filters 
by Eqs. (|8]l to ( 1 1, with a single control parameter w. Figure[T] 
shows cross sections of the filtered characteristic functions for 
two filter widths, w = 1.2 and w = 1.5. The broad and 
narrow curves correspond to the squeezed and antisqueezed 
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Figure 1 : (Color online) Cross sections of filtered characteristic func- 
tions $n of a squeezed state with V x — 0.2, V p — 5.0. The solid 
line shows $n along the squeezed axis with a filter width w = 1.2, 
the dashed line with w = 1.5. The narrow curves belong to the 
unsqueezed axis. 



axes, respectively. Without regularization, $ grows exponen- 
tially in the direction of the squeezed axis, whereas the filtered 
function $$i is square integrable. 

Cross sections of the resulting nonclassicality quasiproba- 
bilities Pq are given in Fig. [2] For both filter widths, they 
clearly display negativities which have their origin solely in 
the nonclassicality of the squeezed state. The larger the width 
of the filter, the more pronounced the negativities become. In 
practice, the filter width is only limited by the experimental 
uncertainties: It should be sufficiently small to keep the statis- 
tical error at a reasonable level. 



Pa(a; w) 




Re(a) 



Figure 2: (Color online) Nonclassicality quasiprobabilities Pq of a 
squeezed state with V x = 0.2, V p = 5.0 along the squeezed axis. 
The solid curve is obtained with a filter width w = 1.2, the dashed 
one with w — 1.5. The negativities unambiguously visualize the 
nonclassicality of the squeezed state. 



We stress that the decisive point in our procedure is that 
the filter has a non-negative Fourier transform. As a conse- 
quence, in our approach, the negativities of the filtered func- 
tion Pq can only be due to the nonclassicality of the quantum 
state under study. In this respect, our method differs from pre- 
vious regularizations of the P function, for example, the one 
by Klauder fT8l , where the used filter has negativities in its 
Fourier transform. Consequently, negativities of the regular- 
ized function Pq do not definitely display the nonclassicality 



of the state. For example, the filtered P function of a coherent 
state is the displaced Fourier transform of the filter itself. It al- 
ways shows negativities, even for the only classical pure state. 
Recently, the P function of a squeezed state has been regular- 
ized by such a filter [ 19 1. The obtained negativities, however, 
cannot be interpreted as the nonclassicality of the state itself. 

In another experiment, the (in this case regular) P function 
of a single-photon added thermal state has been reconstructed 
by using a rectangular filter lfT4l . This requires prior knowl- 
edge about the state to estimate the systematic error caused 
by the regularization, in order to ensure the significance of 
the nonclassical effects. With the methods introduced here, 
the negativities visualize the nonclassicality without any prior 
knowledge about the state. 

The extension of our methods to several radiation modes 
is straightforward. The nonclassicality displayed by the non- 
classicality quasiprobabilities also includes entanglement. To 
directly display entanglement, negativities of entanglement 
quasiprobabilities can be used 1 20 1 . Together with the present 
method this yields powerful tools for characterizing quantum 
systems, which are useful for various applications. 



IV. SUMMARY 

We have introduced necessary and sufficient conditions 
for the nonclassicality of a quantum state which can be di- 
rectly applied in experiments. Universal nonclassicality fil- 
ters and regular nonclassicality quasiprobabilities have been 
introduced which display the nonclassicality of any quantum 
state without prior knowledge of its properties. We have con- 
structed experimentally useful filter functions which only de- 
pend on a single width parameter. The nonclassical proper- 
ties of a squeezed state have been visualized by negativities of 
regular functions, which is impossible with s-parameterized 
quasiprobabilities . 
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Appendix A: General properties of nonclassicality filters 

The approach we are developing shall be applicable to any 
quantum state, also on the basis of experimental data. This 
was already demonstrated in the main text. In addition, neg- 
ativities of regular functions should prove nonclassicality in a 
one-to-one manner. 

Theorem 1 The P function describes a nonclassical state, if 
and only if the regularized function Pq(@; w) shows negativi- 
ties for a sufficiently large, but finite filter width w. 

Proof: Let us assume that a state given by its characteristic 
function $(/?) is nonclassical. Owing to Bochner's theorem, 
this implies that [ 10 1: 
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3N e N, /3i, 
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, (3n e C : 

Y }--«• 



(Al) 
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Let us take a sequence of filters ft w (f3) which satisfies all the 
properties (a)-(c). Then we have for N and /3i, . . . ,/3/y as 



chosen in Eq. ( Al 



lim det 

w— >oo 



- - ^))^=i,...,iv} =Dn<o. 

(A2) 

Since the determinant is a continuous function, there must ex- 
ist a finite wq > such that 



det{($(ft-/%)^o(A-&)) 



} < 0. (A3) 



Hence, we have found a filter £l Wo (ft) such that the filtered 
state <&q(/3;u>o) is nonclassical. Since its Fourier transform 
Pfi(/3; wo) is a regular function but not a probability density, 
it must show negativities to display nonclassicality. 

The other way around, if represents a classical quan- 
tum state, its characteristic function satisfies the conditions 
of Bochner's criterion. Furthermore, each filter fl w ((3) shall 
have a non-negative Fourier transform, and hence it also 
satisfies the conditions of Bochner. Under these assump- 
tions, it is immediately clear that &n(/3;w) also satisfies 
®w(-(3) = ®w(P) and $^(0) = 1. Moreover, the matrix 
($n(A — Pj\ w))i,j=i...Ar is the Hadamard product of the ma- 
trices Pj))i,j=i...N and (O, w (0i-0j))i,j=i...N- If the 
latter two matrices are positive semidefinite, their Hadamard 
product is also positive semidefinite. Consequently, <f>n(/3; w) 
satisfies the conditions of Bochner's theorem. Hence, for any 
classical state with a positive semidefinite P function, the reg- 
ularized function Pq((3; w) is a classical probability distribu- 
tion showing no negativities. ■ 



Appendix B: Decay properties of autocorrelation filters 

Here, we prove the following lemma, which has been used 
for introducing nonclassicality quasiprobabilities. 



Lemma 1 Let u be a real positive number and a real 
function which satisfies C — ||w(/3)e 11 '^' H2 < 00, where ||-|| 2 
is the L 2 -norm. Then the autocorrelation function ofu}(/3), 



w(/3)w(a + /3)d 2 /3, 



(Bl) 



satisfies \\fl(a)e v ' a ' ||a < 00 for any real v < \. 



Proof: The autocorrelation function can be rewritten in the 
following way: 



fi(a) = J uj{P)e um u>(a + f3)e ula+Pl (B2) 

xe -u\2/3+ a \ 2 /2 d 2p e -u\a\ 2 /2^ 

It is bounded from above by 

|Q(a)| < \\co(P)e u ^ 2 u J (a + P)e u ^ 2 e- u ^ +a \ 2 / 2 \\ 1 

(B3) 



xe 



-u\a\ 2 /2 



with || • ||i being the L 1 -norm. Applying Holder's in- 
equality ED in the form < H/llalblbll^lloo with 
f(f3) = u>(t3)e u W\ g((3) = u(a + f3) e u \ a+ ^ and 

Ug-^+al 2 ^!^ = lgiyes 

|n(ot)| < C* 2 e" u|Q|2/2 . 
Since C is finite, we have 

||0(a)e , ' |a|2 ||2 < C 2 ||e^-"/ 2)|a|2 || 2 



(B4) 



(B5) 



where the right-hand side is finite if v < % , ■ 
Remark. If ui((3) satisfies the requirements of lemma[T|for 
all u > 0, the same holds also for fi(/3). 
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